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We present a semiclassical study ol level widths for a class of one-dimensional potentials in the 
presence of an ohmic environment. Employing an expression for the dipole matrix element in terms 
of the Fourier transform of the classical path we obtain the level widths within the Golden rule 
approximation. It is found that for potentials with an asymptotic power-law behavior, which may 
in addition be limited by an infinite waU, the width that an eigenstate of the isolated system acquires 
due to the coupUng to the environment is proportional to its quantum number. 



I. INTRODUCTION 

Semiclassical approaches were essential at the advent 
of quantum mechanics and have since remained a priv- 
ileged tool for learning this subject, for developing our 
physical intuition on new problems and for performing 
analytical calculations [|l|. In particular, in the one- 
dimensional (Id) case the semiclassical (WKB) approxi- 
mation treats a wide class of confining potentials V{q) in 
the limit of large quantum numbers. We can write down 
closed expressions for eigenenergies and eigenfunctions, 
and for power-law potentials the scaling of the eigenener- 
gies with the quantum numbers and the classical actions 
can be established |2||^. Once we have an almost com- 
plete solution of a Id problem we can use it as the basis 
for understanding more complicated physical situations. 
This is the approach of this work, where we analyze the 
problem of a particle in d = 1 coupled to a dissipative 
environment. 

For an isolated quantum system, the ground state and 
the excited states are stable and the corresponding den- 
sity of states is a series of delta functions. Any coupling 
to external degrees of freedom usually renders the ex- 
cited states unstable, thus the levels acquire a finite life- 
time Q which is related to the level width through the 
time-energy uncertainty relation. The density of states 
of the system including the additional degrees of freedom 
still consists of (closely spaced) delta functions, however 
the reduced density of states corresponding to the system 
alone is smeared out compared to that of the uncoupled 
case. 

An example of such a situation is an atom coupled 
to the electromagnetic modes of the radiation field. If 
the atom were isolated from the field there would be no 
transitions from excited states to states lower in energy. 
This changes if we take the coupling to the radiation field 



into account. Then transitions between states may oc- 
cur, and the finite lifetime of the excited states broadens 
the spectral lines associated with the transition. 

Two types of transitions are caused by the coupling to 
the radiation field: stimulated processes and spontaneous 
emission. In the first case photons have to be present, 
while in the second case the mere presence of field modes 
is sufficient to allow for transitions to energetically lower 
states accompanied by the emission of a photon. In the 
limit of zero temperature, the mean number of photons in 
the radiation field vanishes and the broadening of atomic 
levels is entirely due to spontaneous emission. 

In the following, we replace the atom with a particle 
moving in a one-dimensional potential. The environment 
is modeled by a set of harmonic oscillators coupled bilin- 
early to the particle H. Using semiclassical methods, 
we extend the scaling of Refs. H,^ to level widths and 
demonstrate that, for a wide class of confining poten- 
tials, level widths are simply proportional to quantum 
number. To establish this general result we first present 
the theory for the level widths (Sec. |l|) and then apply it 
to a few specific examples: harmonic oscillator (Sec |l|). 
Id box (Sec. |rV), h alf-oscillator (Sec. |Vl[ ) and Coulomb 
potential (Sec. VIl[ ). The semiclassical approximation 
for level widths (Sec. M) simplifies the calculations in 
the examples considered and can be applied to general 
power-law potentials (Sec. |VIII ) . In the concluding chap- 
ter (Sec. IX) we analyze the experimental implications of 
our findings and their possible extensions to higher di- 
mensions. 



II. GENERAL EXPRESSION FOR THE LEVEL 
WIDTHS 

As our model we consider a particle of mass M mov- 
ing in a one-dimensional potential V{q). The spectrum 



of the corresponding Hamiltonian 



i?a = 



2M 



V{q) 



(1) 



is assumed to consist of a discrete part at low energies 
which may be foUowed by a continuous part at higher 
energies. It is on the discrete part of the spectrum (con- 
sisting of eigenenergies En, n = Q, 1,2,. . .) where we fo- 
cus our analysis. To be consistent, we always denote the 
ground state by n = although this may lead to a slightly 
unusual notation. 

Since we consider the limit of large quantum numbers 
we require that the number of discrete eigenstates is in- 
finite or at least can be made sufficiently large. This 
includes for example the radial part of the Coulomb prob- 
lem but excludes the Morse potential. 

The eigenstates acquire a finite width if we weakly cou- 
ple the particle to environmental degrees of freedom. We 
assume that the environment consists of a set of harmonic 
oscillators coupled bilinearly to the particle. This leads 
to the full Hamiltonian 



H = Hs+J2 



j=o 



2711 i 



m^J] 



2*? 



(2) 



implying a coupling between system and environment 
through the Hamiltonian 



j=o 



(3) 



By eliminating the environmental degrees of freedom, we 
obtain an effective operator equation of motion MM , 

/■* 1 dF 1 



with damping kernel 



ry{t) = — / ^^COs(wi) 



M Jq tt oj 
spectral density of bath oscillators 
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'^H-'^E^;;r-^('"-'"^) 
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2mjUjj 



(5) 



(6) 



and an operator- valued fluctuating force ^{t) which we 
do not need to specify further. 

Of great importance is the special case of J(w) ~ M^uj. 
Then the damping kernel becomes ^{t) = 2j6{t). Noting 
that in Eq. (m only half of the delta function contributes, 
the second term becomes jq{t) describing the well-known 
classical damping proportional to the particle velocity. 
This type of damping is often referred to as ohmic be- 
cause such a damping term appears also in equations de- 
scribing electrical circuits containing an ohmic resistor. 



The Hamiltonian (g) provides a microscopic model for 
dissipation in quantum systems in the sense that dissipa- 
tion is due to coupling to additional degrees of freedom. 
It does not, however, pretend that in a real resistor we 
can identify environmental oscillators microscopically. 

Without going into details, we mention that the Hamil- 
tonian (0), besides the fact that the environment can be 
eliminated analytically, provides a good description of 
many realistic systems. It has been well studied over the 
years M and more recently became known as Caldeira- 
Leggett model B in the context of macroscopic quantum 
phenomena. 

Using the Hamiltonian (0) and assuming weak coupling 
between the particle and its environment, we calculate 
the zero temperature width of the n-th level by means of 
the Fermi golden rule 

2 °° 
r„ = — Y, \{m,lj\Hi\n,0,)\^5{En-Era-niUj). (7) 

m,j=0 

This expression describes the decay of the state n to an 
energetically lower state m by one excitation of the j-th 
environmental mode that changes its occupation number 
from to 1. Inserting the dipole matrix element 



'^^l^^l«^>->2™,.. 



1/2 



(8) 



of the j-th environmental oscillator we get 

Tn^n y -^\dnra\''S{En-Ern-hiOj) (9) 

where 



mjLOj 

771, J— 



{m\q\n) 



(10) 



is the dipole matrix element of the system. The proper- 
ties of the environmental modes appearing in Eq. (0) can 
be expressed in terms of their spectral density (0) , so we 
write the level width as 



*. / ^ I ^71777 I '^ 

m=0 



E„ — E„ 



(11) 



This result is valid for arbitrary bath density. Employ- 
ing a cubic frequency dependence for J(w) would lead 
to the natural decay width of an excited atomic state 
due to spontaneous emission, apart from prcfactors aris- 
ing from a proper treatment of the polarization of the 
emitted photons. 

While the following calculations could in principle be 
performed for arbitrary bath densities, we will consider 
only the important case of ohmic damping where the level 
widths are 



2Mj 



77-1 

E 

777 = 



1^77777 1 {En — Em) 



(12) 



The sum over the system eigenstates is restricted since 
an environment at zero temperature cannot excite the 
system into states of higher energy. 

Eq. (|l^ ) constitutes the starting point for a calculation 
of level widths that will be performed in the following 
sections. It represents, up to the factor 7, a sum over 
oscillator strengths 



/, 






' {En — Ejn ) ■ 



(13) 



The finiteness of the upper limit prevents us from eval- 
uating the level widths by the standard Thomas-Reiche- 
Kuhn sum rule [|| for oscillator strengths, 






In 



(14) 



III. HARMONIC OSCILLATOR 

The maybe simplest example is that of a harmonic po- 
tential with frequency ujq. The dipole matrix element in 
this case couples only nearest neighbors, 



2Mujo 



{Vn+ 1 5rn,n+l + \fn 8m,n-\) ■ (15) 



This leads to a simple and well-known result for the level 
widths of a damped harmonic oscillator |lfl| 



r„ = 77,7. 



(16) 



As an illustration we discuss the density of states for a 
damped harmonic oscillator. In view of the relation be- 
tween the partition function of a canonical ensemble, 



m) 



dEp{E) exp(-/3£;) , 



(17) 



and the density of states p{E), one may obtain the latter 
by inverse Laplace transformation pT[ from the partition 
function. For a damped harmonic oscillator, the parti- 
tion function is ||l^ 



Z(/3) = 



1 



n 



^0(^0 ^J-^ vt 



Vnli^n) +t^o 



2 • 



(18) 



Here, 7 denotes the Laplace transform of the damping 
kernel (^) and Vn — 27rn/hp are the Matsubara fre- 
quencies. The product in Eq. (|l8|) does not converge 
in the case of an ohmic environment where jii'n) = 7 
and we therefore introduce a cutoff in the spectral den- 
sity of states. After inverse Laplace transform, according 
to (f^, we arrive at the density of states p{E), for which 
the ohmic limit can be computed if the (infinite) ground 
state energy of system plus environment is subtracted. 
The result is shown in Fig. |l| for ohmic damping with 



j/2ujo = 0.1, omitting the delta function corresponding 
to the stable ground state. The levels broaden with in- 
creasing energy and for large energies only the average 
density of states p = I/Tiloq is seen. 
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FIG. 1. Density of states for a harmonic oscillator 
of frequency uao coupled to an ohmic environment with 
7/26J0 = O.L The delta function corresponding to the sta- 
ble ground state is not shown. 

According to Eq. (|6|) the level widths for the harmonic 
oscillator are proportional to n or, equivalently, to En- In 
general, however, E^ is not proportional to n. Then the 
question arises whether the widths are simply determined 
by n, or by En and whether or not there is a universal 
behavior. 



IV. ONE-DIMENSIONAL BOX 

A more non-trivial example is that of a particle in a 
one-dimensional box of length L. The eigenenergies are 



En 



2Mi2 



(n+l)^ 



(19) 



and dipole matrix elements can be easily calculated. For 
reasons of symmetry a dipole coupling can induce tran- 
sitions only between states of different parity. For n — m 
odd one obtains 



8L (n + l)(m + l) 

7r2 [(?i + 1)2 - (to + 1)2]2 



(20) 



We are interested in the universal behavior for large 
quantum numbers, where d„„i decreases as (n — ?Ti)2. 
Therefore, we define I — n — m and write approximately 
for I <ti m,n 



_ 2L1 



(21) 



We can now evaluate the level widths for large quan- 
tum number n by assuming that the sum over all fi- 
nal states m converges fast enough that only differences 



I <§i n,m are important. Working to leading order in 
l/n we extend the sum over / to infinity and obtain, with 
Eq. (P), 



r -^ 



1 



^^ (2k- 

fc=0 ^ 



1)3 



(22) 



The sum can be expressed in terms of Riemann's zeta 
fmiction as (7/8)C(3) = 1.051 . . . |l3). The level widths 
at large quantum numbers are therefore given by 



r,i,sc = 0.852. . .771. 



(23) 



The subindex "so" stands for "semiclassical" since, as 
seen in the next section, Eq. (^) results directly from a 
WKB approximation. The quality of our approximate re- 
sult is seen in Fig. || where we show the ratio between the 
exact width r„ and the semiclassical width r„.sc- The 
relative error drops below 1% at about n — 38. 

The result (23) indicates that the proportionality in 



the case of the harmonic oscillator should be interpreted 
to be to the quantum number and not to the eigenenergy. 
The prefactor resulting for the Id box is different from 
that of the harmonic oscillator. Therefore a potentially 
universal behavior can only concern the scaling with the 
quantum number, while the prefactor will be specific to 
the given potential. 
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FIG. 2. Ratio between the result ( |23[ ) and the exact width 
r„ for a particle in a box as a function of the state n. 



V. SEMICLASSICAL EXPRESSION FOR THE 
DIPOLE MATRIX ELEMENT 

In the two examples treated above, the harmonic os- 
cillator and the Id box, the wave functions are known 
and the dipole matrix elements are readily calculated. 
To treat generic potentials in the limit of large quantum 
numbers we use a semiclassical approximation for the 
dipole matrix element. Such an ex pre ssion (see Eq. (p3) 
below) was first derived by Debye |l4| and we present it 
in this section to emphasize the physical assumptions on 
which it is based. 



For a Id confining potential V{q) the WKB approxi- 
mation 1 15 1 for the wave function with eigenenergy En 
{n > 1) is 



V'n,sc(g) 



M 



Tp{q,En 



1/2 



1 TT 

COS ( -S{q,En)- - 



(24) 



Here, q lies between the two turning points q\ and 52 
{y(ci\) = ^(92) = £'n, q\ < 92) which determine the lim- 
its of the classical motion at the energy En- The classical 
action is computed from qi to q according to 



Siq,En) 



dqp{q,En) 



with 



p{q,En) = [2M{En-V{q))] 



1/2 



(25) 



(26) 



The quantization condition determining En depends on 
the action along a period of the classical motion 



S{En)^2S{q2,En)^2TTh(n+^ 



and the period T is given by 



T{E) = ^ = 
^ ' dE 



«2 2M 
'^^Q~} ^ 

91 P(q,En) 



(27) 



(28) 



The commutation relation {i/'h)[H,q] = p/M allows 
us to express the dipole matrix element in terms of the 
momentum matrix element. We can then write 



(^n.m. — 



MiEn - Em) 



+ 00 



d 



(29) 



Under the assumption I = n — m -^ n,m, both eigenen- 
ergies have almost the same classically allowed region. 
Neglecting the contribution of the classically forbidden 
region, we restrict the integration to the interval (gi, (72)- 
Moreover, we assume that V{q) is a smoothly varying 
function, and thus the derivative with respect to q is car- 
ried out only in the cosine, which has a small factor Ti 
in its denominator, but not in the prefactor. With these 
assumptions we have 



2h 



T{En — Em) 



«2 



Aq 



-iSiq,En)-Siq,E,n)) 



(30) 



cos 



iSiq,En) + S{q,Em)) 



The second term in the integral is rapidly oscillating and 
is negligible upon integration. For the first term we use 
the relation 



Siq, En) - Siq, E^) ~ tiq, E^){En - E„ 



(31) 



where t{q, En) is the time for a classical particle to move 
from qi to q at energy £"„. With Eqs. (|2^) and (|2^) we 
can write 



dE _ 2Trh 



En — Er^ 

and then find from Eq. (f50| 



1 /"'^ / t(q) 

dn,n-i = 7 dqsin 27r/— ^ 

.T(E) 



(32) 



1 /"'^-^^ / t 

^ j^ dtq{t, E) COS [2^1 TiE) ) ^^^^ 



T{E) 



where the second line has been obtained by means of the 
substitution dq = qdt and partial integration. 

Eq. (|3^) constitutes the interesting result found by De- 
bye which allows to relate the dipole matrix elements in 
the semiclassical limit to the Fourier components of the 
classical motion of the particle. It extends Bohr's corre- 
spondence principle which states that in the limit ?i — > 
the frequencies of an atomic transition should agree with 
electrodynamics. By virtue of Eq. (|33|), a correspondence 
between quantum theory and electrodynamics can not 
only be established for frequencies but also for intensi- 
ties which according to Eq. jl^ ) are related to the square 
of dipole matrix elements. 

Eq. ( p3| ) is the leading order approximation in h. 
Higher order corrections have been derived in Ref. pGJ; 
however we will make no use of them in this work. 

In the presence of hard walls the WKB approximation 
is still applicable and only a change in the phase of the 
semiclassical wavefunction (|2J) is needed to recover the 
result for the Id box. Then the semiclassical approxi- 
mation to the dipole matrix clement is the same as for 
smooth potentials (Eq. (p3)). The period of a classical 
trajectory with energy E in the box is 



TiE) 



/2M 
V E 



1/2 



(34) 



leading in the limit of large n to the forms (|2l] ) and ( [2^ ) 
for the matrix elements and the level widths, respectively. 



VI. HALF-OSCILLATOR 



may then be evaluated by expressing the Hermite polyno- 
mial with the higher quantum number n > m hy means 
of the Rodrigues formula, 



^iivn. — 



1/2 



1 



^MluttJ 2"+"[(2n + l)!(2TO+l)!]i/2 
xj d^eiJ2™+i(e)^|2^exp(-^2)_ (35) 

After repeated partial integration and use of 



— {^H2„.^l{0) 



( ^^,n + l-fc/2 2""^fe(2m+l)! 
^=0 (m+l-A./2)! ' 

(36) 



which holds for k even and yields zero for k odd, one 
arrives at 



^nni — ^ 



MtTU! 

x(-l) 



1/2 



(2TO+1)! 



22(m+n+l)(2n-hl)! 



1/2 



m+1 
n-m+1 \^ 92'/. 



[2(n-0]! 



1=1 



(m - I + iy.{n - ly. 

The sum may be evaluated to obtain 



(37) 



/ n y/^ (_i)n-™+i [(2n + l)!(2m + 1)!]^^^ 



\Mttuj 



On+rn— 1 



nlml 



4(n — ?7i)2 — 1 



(38) 



A corresponding semiclassical evaluation of the dipole 
matrix element using Eq. (p3) leads to 



^n n — / 



fl 

JluJ 



1 

n 

4 



1/2 



1 



4?2 



1 



(39) 



For n ^ / this agrees with the exact result up to an 
irrelevant sign. For the level widths we then obtain 



:;^n ~ 0.810 . . .771 



(40) 



and thus again a proportionality to the quantum number. 
Even though the eigenenergies of the one-dimensional 
box and the half-oscillator depend differently on the 
quantum number and in spite of the fact that the dipole 
matrix elements of oscillator and half-oscillator are quite 
different, the final result for the level widths is propor- 
tional to the quantum number n in both cases. 



Our next example is the harmonic oscillator with an 
infinite wall at g = 0. The particle is therefore confined 
to g > 0. In contrast to the usual harmonic oscillator, 
the dipole matrix element no longer couples only nearest 
neighbor states. The eigenstates of the half-oscillator are 
given by the odd eigenstates of the harmonic oscillator 
with the prefactor adjusted to account for the restricted 
interval of normalization. The dipole matrix element dnm 



VII. RADIAL PART OF COULOMB POTENTIAL 

A less trivial dependence of the dipole matrix elements 
on quantum number is found for our last example, the 
radial part of the three-dimensional Coulomb problem 
for vanishing angular momentum. In this case, an in- 
finite potential wall again restricts the particle to posi- 
tive q where it moves in an attractive Coulomb potential 



—A/q. Unlike the previous examples, the spectrum of 
the Coulomb problem contains a continuous part at pos- 
itive energies. The discrete part of the spectrum consists 
of the eigenenergies 



En 



MA^ 



1 



2%^ {n+lf 



(41) 



Dipole matrix elements can be obtained by solving the 
Schrodinger equation in momentum representation |17| ] 
and read lial 



■/„ 



with 
and 



Sn^ [(n+l)(m+l)]5/2 
M7rA(n + l)2-(m + l)2-'"™ 

■usin[/„(u) - fm{u)] 



Au 



[1 + (n + l)2?i2][l + (TO + l)2y: 



^2„21 



fn{u) = 2(n + l)arctan[(n + l)u]. 
For n,m ^ 1 and n > m this becomes llSt] 



(-1) 



+,2^^^3^^^ /2\ {nmr 



MttA 



1/6 



3/ (n^ — rn2)5/3 



(42) 



(43) 



(44) 



(45) 



It is instructive to derive this result from the semiclas- 
sical expression (|3|). Even though the trajectory q{t) 
cannot be obtained explicitly, we can obtain the matrix 
elements by observing that the behavior of the Fourier co- 
efficients of q{t) is dominated by the most singular part 
of the original function which in our case is the motion 
around the reflection point at q = 0. 

Energy conservation leads to the trajectory at energy 
E in implicit form as 

^^ '^' \E^t +1 = arcsin(,') - q''/\l - q'f' 

(46) 

with q' ~ {\E\/A)q. From this result we compute the 
period 



T^TT 



mV'^ a 



|£;|3/2 



(47) 



and by expanding the arcsin we find, for the trajectory 
around the reflection point encountered at i = 0, 



,= ('J-X"t''^ 



\2M, 
We then write the trajectory for < i < T as 

1/3 / , ^ 2/3 



.(2A\-e"l,_L 



\2M 



T 



'Zrcg 



(48) 



(49) 



where the first term reproduces the most singular behav- 
ior of the trajectory around i = and t — T while the 
second term contains the more regular part and is of no 
further interest. From Eq. ( [f9|) we find with Eq. (|3^ ) for 
ri S> / 3> 1 and up to an irrelevant sign 



31/6 



MnA 21/3 



n 



(50) 



Written in symmetric form this agrees with Eq. (|45| ) in 
the limit n ^ I. 

The level widths can now be obtained from Eq. (|lj) 
with (|l]) and (|^) as 

6V3[r(2/3)]% /7\ 
r„,sc = 2 C o 7" = 0-477 . . . 7n . (51) 



Again, we find the proportionality to the quantum num- 
ber n which, in view of the complicated structure of the 
dipole matrix element, may come as a surprise. This 
indicates that this proportionality might be more uni- 
versal. We will indeed prove it in the following section 
for power-law potentials with and without walls. 



VIII. SEMICLASSICAL TREATMENT FOR 
POWER-LAW POTENTIALS 

In this section we generalize the previous results to 
power-law potentials of the form V{q) = Aj^j". The am- 
plitude A and the exponent a should have opposite signs. 
The possible coordinates might be limited by an infinite 
potential wall, and we assume that such a wall is present 
at g = whenever a < 0. The case a = will be ex- 
cluded because it requires two walls and thus becomes 
the box potential already treated in Sec. 



IV 



Further- 
more, we have to restrict the exponent to a > —2. At 
a — —2 the action becomes independent of energy as 
will become clear from Eq. ( p6| ) below so that we have 
to exclude this pathological case. We emphasize, that an 
attractive 1/g^ potential can never appear in an radial 
equation of motion in d > 1 where the angular degrees 
of freedom have been eliminated. 

Since our semiclassical approach requires sufficiently 
large energies, the discussion applies also to potentials 
which effectively behave like a power-law potential at 
higher energies, e.g. the quartic double-well potential. 

According to Eq. (|lj) the n dependence of the level 
width is determined from the dipole matrix element dnm 
and the energy difference {En — Em)- Within the hy- 
pothesis established in Sec. M (which was more general 
than for power-law potentials) we saw that (£"„ — E„i) is 
simply inverse proportional to the period of the classical 
motion between the turning points (Eq. (p2|)). 

To evaluate the dependence on n of the level widths 
(n2) we employ scaling properties of the energy conser- 
vation condition 



Introducing a diniensionless coordinate 



q' = 



E 



l/c 



as in ( [46[ ) and diniensionless time 



(p2) becomes 



f' — f 

Ml/2 



1 = -g'2 + sign(A) q"^ 



(52) 



(53) 



(54) 



(55) 



where q' = dq' /At' and signM) denotes the sign of A. 
The quantization condition (|27|) reads in scaled variables 



S{E) = 



j^jl/2^{2+a)/2a 
\A\^ 



dq' q' ~ 27r?in , 



(56) 



where the integral runs over one period. On the right- 
hand-side we have omitted an n-independent term which 
depends on the presence of walls in the potential. With 
S' = f dq'q' we then find for the energy eigenvalues 



E = 






\A\_ 
M" 



l/(2+a) 



2q/(2+q) 



(57) 



This is in agreement with previous results for power-law 
potentials |g,^ where S' has been evaluated explicitly. 
On the other hand, Eq. (57) is still correct for sufficiently 
large n if the potential behaves only asymptotically like 
a power-law. Then, in general, S' can no longer be eval- 
uated analytically. 

The semiclassical dipole matrix element (|3|) reads af- 
ter scaling 



Ojn. .rt. — / 



E 



1/q 



ra' '■ 



with 



d; = -4 



t'iq') 



TTlJ^, \ T' 



(58) 



(59) 



where T' denotes the period T scaled according to (M| 
In view of Eqs. (H), (f^) and (||), the level width^ 
can be expressed as 



8n^2 + a 
S'^ 2a 



j:idf 



(60) 



Apart from the factor n this result still depends on the 
state number n via the upper limit of the sum. As a 



last step, we therefore have to consider the convergence 
properties of this sum. 

For a > and in the absence of an infinite potential 
wall, the particle velocity is always continuous and, as a 
consequence, the scaled dipol moment (p9[) will decay at 
least as l^^. This still holds for a wall with finite poten- 
tial to its one side because in the worst case the refiection 
will lead to a triangular cusp singularity in the trajectory. 
In this case, the dipol moment will decay as l~^. 

The case of negative exponent a with an infinite po- 
tential wall at the divergence at g = is more interesting. 
Close to q = wc may neglect the constant on the left- 
hand-side of (|55|). Assuming that the reflection happens 
at t' = we find for the trajectory close to the reflection 
point q' ^ |i'p/'2~"). Proceeding as in Sec. VII, one finds 
for the scaled dipole matrix element dj ~ ;-(4-a)/(2-a)^ 
For a = — 1 this is in agreement with (pO|). For a > —2 
it follows that dj decays always faster than l^^^. 

As a consequence, the argument of the sum in Eq. ( |60| ) 
decays faster than 1/P for all potentials under considera- 
tion. Neglecting terms of order 1/n, as is consistent with 
our previous approximations, we may extend the upper 
summation limit to infinity and arrive at our final result 
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(61) 



For sufficiently large energies, the level width is therefore 
proportional to the state number n. We point out that 
the proportionality constant depends on a and 7 but not 
on M and A. While these properties are special to the 
case of ohmic damping, an extension to other bath den- 
sities would be straightforward along the lines presented 
here. 



IX. CONCLUSIONS 

Semiclassical methods for radiative lifetimes, ubiqui- 
tous in spectroscopic measurements, have mainly been 
derived for excited states of hydrogenlike atoms pi[ . 
Here, we applied such methods to a whole class of one- 
dimensional potentials which behave like a power-law for 
large energies and may contain one or two walls. Using 
a semiclassical expression for the dipole matrix element 
in terms of the Fourier transform of the classical paths 
we have shown that for ohmic damping the level width is 
proportional to the number n of the state. This result is 
valid for sufficiently large n, therefore requiring weak cou- 
pling to the environment in order to be observable. The 
proportionality of level widths with the state number of 
the harmonic oscillator could thereby be generalized to a 
large class of one-dimensional potentials. The prefactor 
of the linear law depends on the specific potential. 

Extensions of our results to higher dimensions are in- 
teresting since we will be approaching experimental re- 
ality and will be able to investigate the relevance of the 
integrability of the classical system on level widths. The 



chaotic or integrable character of the classical dynam- 
ics has been shown to be of crucial importance for other 
observables like the orbital magnetic response [ pO[ . 

The problem of an integrable 2d system, where we can 
separate the dynamics, can be reduced to two Id prob- 
lems. For the 2d box it is readily demonstrated that the 
mean behavior of r„ is again a linear increase with n. 

The problem becomes more involved in chaotic systems 
where the energy is the only conserved quantity and the 
explicit forms of individual wavefunctions are not known. 
When increasing the number of degrees of freedom of the 
system from 1 to 2 we start seeing fluctuations of level 
widths with respect to a secular behavior. Since level 
widths are given by matrix elements of the dipole oper- 
ator we expect that the fluctuations will be more impor- 
tant in the integrable case, where selection rules impose 
wide variations according to the quantum numbers of the 
initial and final states. 



ACKNOWLEDGMENTS 

We have benefited from discussions with J.-Y. Bigot, 
H. Grabert, S. Kohler, and S. Otto. Part of this work 
has been carried out while one of us (GLI) was at the 
Centre d'Etudes de Saclay with financial support from 
the Volkswagen-Stiftung. 



[1] Matthias Brack and Rajat K. Bhaduri, Semiclassical 
Physics, Frontiers in Physics Vol. 96 (Addison- Wesley, 
Reading, 1997). 

[2] Uday P. Sukhatme, "WKB Energy Levels for a Class of 
One-Dimensional Potentials," Am. J. Phys. 41, 1015- 
1016 (1973). 

[3] J. F. Cariiiena, C. Farina, and Cassio Sigaud, "Scale in- 
variance and the Bohr- Wilson- Sommerfeld (BWS) quan- 
tization for power law one-dimensional potential wells," 
Am. J. Phys. 61, 712-717 (1993). 

[4] V. Weisskopf and E. Wigner, "Berechnung der 
nattirlichen Linienbreite auf Grund der Diracschen Licht- 
theorie," Z. Phys. 63, 54-73 (1930). 

[5] see e.g. Thomas Dittrich, Peter Hanggi, Gert-Ludwig In- 
gold, Bernhard Kramer, Gerd Schon, and Wilhelm Zw- 
erger. Quantum transport and dissipation (Wiley- VCH, 
Weinheim, 1998), Chap. 4. 

[6] V. B. Magalinskii, "Dynamical Model in the Theory of 
the Brownian Motion," Zh. Eksp. Teor. Fiz. 36, 1942- 
1944 (1959) [Sov. Phys. JETP 9, 1381-1382 (1959)]. 

[7] see e.g. Robert J. Rubin, "Statistical Dynamics of Sim- 
ple Cubic Lattices. Model for the Study of Brownian Mo- 
tion," J. Math. Phys. 1, 309-318 (1960) and 2, 373-386 
(1961); I. R. Senitzky, "Dissipation in Quantum Mechan- 
ics. The Harmonic Oscillator," Phys. Rev. 119, 670-679 
(1960) and 124, 642-648 (1961); G. W. Ford, M. Kac, 



and P. Mazur, "Statistical Mechanics of Assemblies of 
Coupled Oscillators," J. Math. Phys. 6, 504-515 (1965); 
P. Ullersma, "An exactly solvable model for Brownian 
motion," Physica 32, 27-55, 56-73, 74-89, 90-96 (1966); 
Robert Zwanzig, "Nonlinear Generalized Langevin Equa- 
tions," J. Stat. Phys. 9, 215-220 (1973). 
[8] A. O. Caldeira and A. J. Leggett, "Quantum Tunneling 
in a Dissipative System," Ann. Phys. (N.Y.) 149, 374- 
456 (1983). 
[9] W. Thomas, "Uber die Zahl der Dispersionselektronen, 
die einem stationaren Zustande zugeordnet sind," Natur- 
wissenschaften 13, 627 (1925); W. Kuhn, "Uber die 
Gesamtstarke der von einem Zustande ausgehenden Ab- 
sorptionslinien," Z. Phys. 33, 408-412 (1925); F. Reiche 
and W. Thomas, "Uber die Zahl der Dispersionselektro- 
nen, die einem stationaren Zustand zugeordnet sind," 
Z. Phys. 34, 510-525 (1925); for a presentation in En- 
glish language see e. g. C. Cohen- Tannoudji, B. Diu, and 
F. Laloe, Quantum Mechanics (Wiley, New York, 1978), 
Vol. 2, p. 1318. 

[10] see e. g. Claude Cohen-Tannoudji, Jacques Dupont-Roc, 
and Gilbert Grynberg, Atom-Photon Interactions (Wiley, 
New York, 1992), p. 326. 

[11] A. Hanke and W. Zwerger, "Density of states of a 
damped quantum oscillator," Phys. Rev. E 52, 6875- 
6878 (1995). 

[12] see e. g. Ulrich Weiss, Quantum Dissipative Systems 
(World Scientific, Singapore, 1993), p. 51. 

[13] Handbook of Mathematical Functions, edited by Milton 
Abramowitz and Irene A. Stegun (Dover, New York, 
1972), p. 807. 

[14] P. Debye, "Wellenmechanik und Korrespondenzprinzip," 
Physik. Zeitschr. 28, 170-174 (1927). 

[15] Martin C. Gutzwiller, Chaos in Classical and Quantum 
Mechanics (Springer, Berlin, 1990). 

[16] Robert Karrlein and Hermann Grabert, "Semiclassical 
theory of vibrational energy relaxation," J. Chem. Phys. 
108, 4972-4983 (1998). 

[17] Eugene V. Ivash, "A Momentum Representation Treat- 
ment of the Hydrogen Atom Problem," Am. J. Phys. 40, 
1095-1097 (1972). 

[18] S. M. Susskind and R. V. Jensen, "Numerical calculations 
of the ionization of one- dimensional hydrogen atoms us- 
ing hydrogenic and Sturmian basis functions," Phys. Rev. 
A 38, 711-728 (1988). 

[19] S. P. Goreslavskii, N. B. Delone, and V. P. Krainov, 
"Probabilities of radiative transitions between highly ex- 
cited atomic states," Zh. Eksp. Teor. Fiz. 82, 1789-1797 
(1982) [Sov. Phys. JETP 55, 1032-1036 (1982)]. 

[20] Klaus Richter, Denis UUmo, and Rodolfo A. Jalabert, 
"Orbital magnetism in the ballistic regime: geometrical 
effects," Phys. Rep. 276, 1-84 (1996). 

[21] H. A. Kramers, "On the Theory of X-Ray Absorption 
and of the Continuous X-Ray Spectrum", Philos. Mag. 
46, 836-871 (1923); K. Omidvar, "Semiclassical formula 
for the radiative mean lifetime of the excited state of the 
hydrogenlike atom," Phys. Rev. A 26, 3053-3061 (1982); 
Hermann Marxer and Larry Spruch, "Semiclassical es- 
timation of the radiative mean lifetimes of hydrogenlike 
states," Phys. Rev. A 43, 1268-1274 (1991). 



